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FOREWORD 
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publication are in preparation and will be released subsequently. 
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INTRODUCTORY REMARKS 


In a variety of situations of engineering interest, materials 
behave in a nonlinear manner, that is to say, the response to excitation 
deviates from a simple linear relationship. In years past, it was 
sufficient to ignore this behavior or, in some instances, to approximate 
it in what is now viewed as a crude manner. More recently, however, 
there has been impetus to look closely at nonlinearities . The motivation 
derives from several sources. On the one hand, design requirements 
have grown more stringent as the result of the need for improved 
performance of materials, the use of more expensive materials, and the 
desire for better prediction of response to high excitation. On the 
other hand, our capabilities have grown not only with the overall 
advances of technology but also and perhaps more importantly, with the 
availability of large high-speed computers. These tools have provided 
researchers and engineers both with the capacity to solve problems 
that once were posed largely for intellectual interest. Computers have 
also served to reduce the disparity of language between researcher 
and practicing engineer; the one can translate complex mathematical 
results into computer programs, and the other can use the code to solve 
problems directly. 

One of the kinds of nonlinearity that is of extensive interest 
occurs in metals and their alloys where excitation is sufficient to 
induce plastic flow. Such behavior is expected in structures designed 
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for high strength and is sometimes observed in structures otherwise 
intended to be quite stiff. Plastic flow also occurs in regions of 
high stress concentrations such as holes, cutouts, and notches, 
sometimes quite unintentionally. The issue in such cases becomes 
that of determining the integrity of the structure after plastic 
flow has occurred. In addition, there is frequent interest in the 
mechanics of the flow process itself partly to determine whether the 
process can either be avoided or be turned to advantage. 

The theory or, more accurately, theories of plasticity are hardly 
new to students and practioners of mechanics. The literature dates 
from the last century, and elements of plasticity are to be found in 
engineering curricula, both undergraduate and graduate, throughout the 
country. Since 1950, the publication date of Hill's book [1]*, 
there has been a resurgence of interest and effort in plasticity. 

Yet, in certain respects, such attention has been accompanied by 
inherent limitations which may be of consequence in a variety of 
applications. The use, for example, of 'perfect' plasticity fails to 
model many materials of technological importance; rigid plasticity is 
inappropriate to structural situations in which elastic and plastic 
strains are both of interest; and slipline techniques tell little of 
the process that precedes and engenders flow. 

Of more recent origin is what we term here the theory of elasto- 
plastic flow. Developed by the author, his colleagues and students 

*Numbers in brackets denote entries in the list of references following 
the text. 
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over the past few years, elasto-plasticity results from a special 
formulation of information already in hand. This formulation is 
believed to serve two purposes. First, it provides means for solving 
a broad variety of problems involving' plastic flow, a few of which are 
indicated in later sections of this report. Second, it suggests an 
approach for formulating theories of other non-linear phenomena. While 
we do not discuss the latter feature in this report, we fully expect 
it to offer a significant potential to other researchers and engineers. 

This report is divided into several sections which give general 
and special forms of the constitutive relations, together with other 
relations which provide the full theory; reductions to two spatial 
dimensions; solutions to elementary problems; numerical methods for 
solving more complex problems; and a review of some solutions that 


have been obtained. 
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CONSTITUTIVE RELATIONS IN GENERAL 


As suggested above, there are several theories of plasticity, 
distinguished by a variety of relationships between strain or other 
kinematic measures, and stress. Accompanying this range, there are 
many approaches available for deriving or developing these relationships, 
and the interested reader may care to peruse some of the texts now 
available to see these approaches. Since our discussion is not as 
heuristic as that found in a text, we proceed on a more axiomatic 
basis, relying on just two primary assumptions. 

In addition, we omit dynamic effects so that the flow process to 
be described is taken to be slow, in the sense that no wave phenomena 
occur. For the most part we deal with infinitesimal strains, although, 
attention is given below to finite strain effects. While body forces 
will appear in the formulation of the next section, they are regarded 
largely in a formalistic sense and are excluded in specific problems. 


The first primary assumption is that there is a loading function 


f which depends upon the stress field* a^, 
and the plastic strain energy density W^. 


the plastic strains 
At each point of the body. 


if f < 0, no flow occurs; where f = 0, flow may occur and we write this 


*Standard indicial notation and its associated conventions are used 

throughout, except where noted. The range of roman indices is 1,2,3 

and 6.. is tne Kronecker delta. 

13 
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condition on the loading function as 



J p ]w (p) ) 

il 


< 0 


( 1 ) 


This statement may be visualized somewhat more easily if we write 
the particular form 

fCa.^ J p) , W Cp) ) = <f>|a i;j ) - W (P) ) < 0 (2) 

which is used throughout most of our development. In (2) ) is 

a function which may be regarded as having two roles. In the hyperspace* 

whose axes are the stresses a... <b is a measure of the vector from the 

il 

origin to the point describing the stress field at a point in the body. 

If this measure is less than the scalar ip, no flow ensues. Altema- 

I 1 

tively, the relation i 

4>O i; j) ~ 4 > = 0 (3) 

is a surface in stress space; if the stress vector is such that its tip 
lies on the surface, flow may occur depending upon how the stress vector 
changes as a result of additional excitation. Note that flow will 
alter ip so that the stress vector will never pierce the surface. 

The second primary assumption is that hypothesized by Drucker [2] . 
Without reviewing the arguments in full detail, we note that the 
hypothesis concerns a prestressed body subject to a small increment in 
load. This self-equilibrated loading is applied and then removed. 


*Altematively, one may think in terms of a 3-space whose axes are the 
principal stresses oi, o 2 , 03 . While the alternate is more easily 
visualized, the form in ( 2 ) is more useful in formulation of the 
constitutive equations . 
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The resultant stress increment Sa„ engenders an increment in plastic 
strains 6e!^ and, largely on thermodynamic grounds, Drucker requires 


that 


So. .6e^ 


0 


il il 

That the quantity on the left of (4) is nonnegative requires that no 
work may be extracted from the process. We may also interpret (4) 
to imply that the motion in the process corresponds to the loading, 
and vice versa. 

The utility of (4) lies in the conclusions that may be drawn 
from it. The first is that the surface described in (3) is convex, 
and the second that the direction of is normal to that surface. 


Hence we have 


5e^ = A3f/3o. . 
il il 


( 4 ) 


(5a) 


or 


= \d<p/do. . 
il il 

depending on whether (1) or (2) is the appropriate form of the loading 
function. The third conclusion is that the increments of stress and 
plastic strain are proportional to one another and, as that result is 
central to our formulation, we consider it in further detail below. 

The arguments necessary to develop convexity and normality may be 
found in Drucker's paper [2] or Fung's book [3], chapter 6. 


(5b) 
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Consider a point in a body whose loading history is such that 
(3) applies. Let an increment in excitation be imposed such that 
flow occurs; it follows that <j> and \p will change in concert: 

6<J> - S\p = 0 

or, to first order terms, 

(3<f>/3o. .)6 ct. . - [3^/3e.^ + (3iJ;/3W ( ' P) )a . . ] 6e = 0 

i y 13 1 3 13 13 

where we have used the differential definition of plastic strain 
energy density in the increment: 

= a. 

1 3 13 

Substituting. (5b) into (6) and solving for X gives* 


X = 


(3^/3a i j)6a i _. 




mn 


mn 


(7) 


so that there results 


13 


(P) _ 


(3<J>/3a i ^) C34>/<5a k ^) 6a k ^. 


+ (9^/3W Cp) )a m J(3<i>/3a m J 


(8a) 


mn 


mn J 


mn 


and the proportionality between plastic strain and stress increments 
is clear. Had we used (5a) , we would find 


13 


(3f/3a.,)(3f/3a k£ )6a k£ 


+ Of/sw'Ph^lOf/aa™,) 


(8b) 


mn 


mn 


*An obvious substitution of dummy indices is required in this step. 
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Either of (8a, b) is often termed the flow rule for incremental plasticity. 

It is useful next to introduce an intrepretation of the function <j). 
As we have described it, it is suggestive of an equivalent stress, and 
we now specify that 

<b = x 

eq 

and it should be noted that eq are not indices. There is then a 
corresponding equivalent plastic strain increment defined by 

x Se^ = a. = 6W^ (9) 

eq eq i] i] 


and we infer a functional relationship between the two equivalent 
quantities. In particular, we have at any value of x such that 
flow is in progress 


(P) 


6 t /6e 
eq eq 


= 2y 


(p) 

eq 


(10) 


It may be observed that, once suitable experimentation is performed 
to find the form of the function <f>(0.jj)> the (scalar) equivalent 
stress and plastic strain relationship may be identified, and that 
the modulus of that relationship may be determined. If we combine 
(8a), (9), and (10) we find that 


(P) 
6e . . 
ij 


K9<»/8a i; j)0<f>/9a k£ ) 

2y ( ' P ' ) (3$/3o )cr 
eq mn mn 


60 


k l 


( 11 ) 
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We have thus a form of the flow rule that depends on <J> and the more 

physical y^ , rather than \p. It does not appear possible to reduce 
eq 

(8b) to an analogous form since f is not dependent only on a ^ . 

The total strain increment comprises both elastic and plastic 
portions, and we have so far dealt with the latter. The former 
follows Hooke's law written in incremental form, viz: 

2 y6e^ = 6a . . - [v/ (1+v) ] 6a. . 6 . . (12) 

ij i] L ^ J 1 kk 13 v ' 

where y is the elastic shear modulus and v is Poisson's ratio. Combining 

(11) and (12) we have the total flow rule 

2y6e„ = 6a i;j -[v/(l+v)]6a kk 6 i;j + r(3<j>/3a i; .) (3<f>/ 90 ^) 60 ^ (13) 

where we have set 

T = [y/y ^ ] <t>/[(3<j>/3a ) a ] 

1 eq J L v mn mn 

as a matter of convenience. The inverse to (13) is 


6a. ,/2y = 6e. . + [v/(l-2v) ]6e. . 6 . . 
ij ij 1 J J kk ij 


( 14 ) 


{34>/3a ij -^[v/(l-2v)]3^/3a rs 6 rs 6. .}{3(j i /3a k ^[v/(l-2v)]3(i)/3a pq 6 pq 6 k£ }6 £u 

(3<}>/3a ){3<j>/3a + [v/ (l- 2 v) ] 3<j>/3a 4 . 6^ 6 } +1/T 

mn Y mn L K tu tu mn 

Terms of the form ( 34>/ ^ a rs ) ^ rs are to be observed in (14). The degree 
to which such terms do not vanish is a measure of compressibility of the 
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plastic strains. That is, in (11), we have the plastic dilatation 
which so far is not presumed to vanish. If it is required, however, 
that plastic strains are incompressible, (14) reduces to 


Sa . ./2y = 6e. .+ [v/(l-2v) ]6e. . 6 . . 
ij ij L J kk ij 


0<j)/3o. j)(3<fr/3o u 


(3<f)/aa )(3d>/3a 

mn 


mn 


)6 V 

) + 1/r 


a somewhat more compact form and one which corresponds to the 
usual assumption (in metals) of incompressibility of plastic strains. 

While (8b) does not reduce to a form analogous to (11), it 
may be inverted to a form similar to (14); the result is 


<5a. ./2y = <Se . .+ [v/l-2v) ] <5e. . <5 . . 
ij ij kk ij 


{(3f/3a. ■) + [v/(l-2v)]3f/3a rs 6 rs 6..H3f/3a k£+ [v/(l-2v)]3f/3a p q3pq6 k£ }6 gk£ 


(3f/3a ) { 3f/3a + [v/(l-2v) ] 3f/3a^ 6^ 6 } +l/r* 

v mn^ mn 1 J tu tu mn 


where 


T* = -2y/{[3f/3e + (3f/3W^^ )a ](3f/3a )} 

L mn mn J mn 


Note that utility of the last expression awaits a definite form for f, 
although one may use this result in a formal sense in the development 
below. 


In each of these cases, we find relationships that show the pro- 
portionality between the increments in stress and strain. Symbolically 
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(15) 

(16) 

where Ch and are a generalized compliance and modulus whose 

specific forms depend upon the loading function selected or otherwise 
taken to model the material. To the extent that other models may be 
put into this same form, we may proceed to develop the complete set 
of equations for elasto-plastic flow, as shown in the next section. 


we may therefore write 


and 


2)jSE ij = C ijk£ S V 


ia ii l2v * hjkZ U kl 
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GOVERNING EQUATIONS 

The derivation of a flow equation of the form of (15) and the 
existence of its inverse (16) permits development of the full theory 
to proceed. We take as a starting point the fact that the constitutive 
equations relate increments of stress and strain in an explicit manner. 

It follows that the remaining elements of the theory should deal with 
incremental quantities directly; to do otherwise leads perforce to 
nonlinearities in the theory and therefore difficulty in solving 
problems. 

Experience has shown that the simplest formulation will be in 
terms of displacement increments rather than potential functions 
(or. their increments), such as any of several types of stress functions. 
Accordingly, we write the strain-displacement equations in the form 

6e. . = (6u. . + 6u. .) / 2 (17) 

il i,l l,i 

where 6u^ are the displacement increments. These quantities will 
become the dependent variables of prime interest. Note that (17) is 
tensorially correct when the differentiation is interpreted as 
covariant inasmuch as the development here presumes small deformations 
(or displacement gradients). This feature characterizes the remaining 
development except where noted in the following sections. 

The equilibrium equations are also written in incremental form. 

6a. . . + 6X. =0 

11,1 * 


Formally, we have 
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where <$X^ are increments in the body force. For most problems of interest, 
however, we drop reference to body force and write simply 


So . . . = 0 

ij >3 


(18) 


as the equations of equilibrium. 

Note that both (17) and (18) are in effect differential equations 
in both space and some time- like parameter that reflects the sequential 
nature of the loading and the process of elasto-plastic flow. It 
follows that (17) and (18) must be accompanied by initial conditions. 

One might specify, for example, the position of the body at an initial 
'time', along with its then-existing state of stress and strain. While 
in many problems, these field quantities will be null, it is clear 
that problems involving, say, a residual stress state may come under' 
consideration; the mathematical need for initial conditions provides 
the proper means for inserting such physical information into the problem. 

We now assemble the foregoing. (17) is inserted into (16) and, 
owing to the symmetry of the generalized modulus, we have 

6 °i) /2v m E ijk£ 5 \,£ fl9 > 

Equilibrium is enforced to obtain 


(E ijM 6 u k,£ ) ,j " 0 


( 20 ) 


which is an extended or elasto-plastic form of Navier's equations. 
It is this (set of three) equation (s) which governs the flow of an 
elasto-plastic medium. 
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Boundary Conditions: It should be reasonably clear that (20) is 

an analogue to the equations of classical elasticity; the differences 
are that the modulus is variable rather than constant, and that the 
variables are the increments of displacement rather than their 
instantaneous values. We may use this analogue to write boundary 
conditions. At the boundary of the domain in which (20) is to be 
integrated, one specifies either increments in displacement <5u^ or 
increments in traction 6t^, the latter being given by an adaptation 
of Cauchy's equations 

st i = 5 o ij n j ■ 12 

having employed (19). In (21) m is the unit outward normal to the 
boundary. It is further possible to specify a mixed boundary-value 
problem in exactly the same manner as is done in elasticity. 

Integral Theorems: Because the theory is formulated in a manner 

to identify incremental quantities as dependent variables, it is 
straightforward to devise integral theorems that may be of use in solving 
problems. Some of these have been described elsewhere, and we repeat 
them for completeness. In order to avoid complexity of notation, we 
replace the incremental notation by that usually associated with rates; 


Su. -*■ u. 
l i 


Sc. 

1 3 i J 

So. . -*■ o. . 
ij 


thus 
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so that the variables of interest are instantaneous rates of change 
(with respect to the inferred time-like parameter noted above) 
rather than increments . 

We have first an analogue to the theorem of minimum potential 
energy; consider the functional 

II = y / E. D e. . e, ^dV - / X.u.dV - / t.u.dS 

y i]kl ij kl lii J 11 

V VS 

a 

in which V is the interior of the domain of interest, S is its 
surface or boundary, and S^ is that portion on which the traction 
rates t^ are specified. It may be shown [4] that, if all displacement 
rates that satisfy stated boundary conditions on S , where 
S = + S u> those which satisfy the equations of equilibrium in rate 

form are distinguished by a stationary value of II; moreover, that 
value of II is a minimum. The procedure by which this theorem is 
established is given in [4] and by the pattern set in [3] , chapter 
10. Basically, one sets the first variation of H to zero, i.e.. 


an 


2 y / E. D z. .<5e. „dV - / X.du.dV - / t.6u.dS 

/ijk£ijk£ lii Jii 


V 


= 0 

and integrates terms in the first integral by parts, using (17). The 
results are (20), modified to include body forces, and (21), in rate 


form. 
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We consider next the functional 


n* = 


l 

4y 


/ 


C ijk£ a ij a k£ dV 


/ ii.X.dV - / ix.t.dS 

11 J ii 


V VS 

u 

Using the same type of procedure, it may be shown that, of all stress 
rates that satisfy the equations of equilibrium rate and boundary 
conditions on S , those which lead to a compatible strain rate field 
are distinguished by a stationary value of n*; moreover, that value 
of n* is a minimum. This theorem is less useful than the one above, 
since its proof and use require introduction of stress (rate) functions 
which we view, as noted above, as awkward. 

Third, we write the functional 


r, f.... 

= J - x i u i) dv - J V E i;f tu i,V u :i,i )/21 


dv 


V 


/ t.u.dS - I n .o . . (u. -u. * 

ii J J !J 1 i 


)dS 


u 


• * 


in which u^ are the displacement rates specified on S^. We compute 

• • • • 

the first variation of P, allowing u^, and to vary in V, u^ 

• • • • ^ 

to vary on S , and o. . to vary on S . Note that X. , t. , and u. are 
7 o lj 1 u 111 

prescribed. Computing 6P, integrating by parts those terms involving 
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> 


and setting the result to zero, we find 


CT ij = 2yE ijk£ e k£ 

in 

V 

a.. . + X. = 0 
1 1,1 i 

in 

V 

e. . = (u. .+u. .)/2 

il i,l l,i 

in 

V 

• • 

t. = cr. .n. 

i il 1 

on 

S a 

• • * 

u. = u. 

on 

S 


11 u 

Thus this analogue to Reissner's theorem gives a recapitulation of the 
various equations derived above. A complementary form of the theorem 
derives from the functional 


p* = 


- At- c. . 

J 4y ij 


jk£ 0 ljV fXu i ldV 


f l '°i 


(u. 
11 1 


+ u . ) / 2 ] dV 

• J j , 


V 


— / u.t.dS - / (u. -u*) a . .n . dS 

J 1 1 ; 1 i i] j 

s s „ 

a u 

but is of limited utility as noted above for n*. 

Finally, we note the existence of a form of Betti's reciprocal 
work theorem. It has been shown [5] that 


/ 

V-V 


a.'.E^dV 
il il 



r.a. .dV 
il il 


V-V 


The primed stress and strain fields are those associated with an 
elastic field; in particular we consider Kelvin's problem of a point 
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load in an infinite elastic body. V' is a ball surrounding the 
load point, deleted from the integrals owing to the singular nature 
of the local stress and strain fields. Note that rates or equivalently, 
increments, form the basis of reciprocity. Due to the inferred 
ellipticity of the governing equations [5,6], the u^ are continuous 
and should possess continuous second derivatives . Thus the divergence 
theorem holds and the reciprocal theorem may be written as 


/ 

S+S" 


t . u. dS 
1 l 



u . t . dS + 

l l 



or.^av 

1J 11 


where is the bounding surface of V". Taking the limit as V' -*■ 0, 
we have 

u. + f T. .u.dS = f\J. .t.dS + (V. e^dV 
l J ij j J i] j ■ J^U k jk 

S S V 


as an extended form of Somigliana's identity. In this relationship, 
note that the plastic strain rate is proportional to the stress rate 
and therefore the total strain rate, through either of (8), and (10). 
Thus this integral statement may be written to relate surface tractions 
and displacements may be eliminated [5]. The kemals T^ . , , and 

^ijk are §^ ven [5] and need not be repeated here. 

Undoubtedly, other relationships analogous to those familiar in 
elasticity may be developed for elasto-plastic flow. The reason for 
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this, and the central feature of the present theory, is the attention 
given to increments (or equivalently, rates) of field quantities as 
the dependent variables of interest when a problem is to be solved. 

What this type of formulation does is to shift the non linearities of 
the physical process into a slightly subordinate position. Mathemati- 
cally, the governing equations are not nonlinear . Instead, they are 
quasilinear; the nonlinear aspects of the process show up as coefficients 
in the governing equations and, as described below, the solution of 
problems is greatly facilitated. Moreover, the equations themselves 
are accessible to study; we have, for example, looked in varying 
degrees of detail at the character of the equations [5,6] and found 
them to be elliptic so long as there is some work-hardening. This 
nature of the equations has enormous implications for the means to be 
used in solving problems. We know, for example, that so long as 
there is some work -hardening, no slip lines will occur. Thus many 
of the solution techniques developed for rigid and perfectly plastic 
materials are demonstrably inappropriate. In their place, however, 
one is free to choose methods appropriate to elastic problems; the 
chief numerical methods thus available are finite elements, finite 
differences, and boundary integral equations. Each of these offers 
certain advantages and disadvantages to the analyst, depending upon 
the particular problem at hand, and he has access to the extensive 
experience relating to each as he proceeds to attack his problem. We 
discuss some aspects of solution strategy in a subsequent section of 
this report. 
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SPECIAL CONSTITUTIVE RELATIONS 

In the penultimate section, we discussed constitutive relations 
in general, arriving at relations of the form 


2l,6£ ij = C ijk£ 6o W 

and the inverse 


(15) 


S \j /2v = E ijW 5E k£ 


(16) 


It is important that relations of this type be used so that the 
remaining development of the governing equations may be followed. 

Within these limits, however, a considerable variety of specific 
constitutive relations may be employed and, in this section, we 
indicate a few that may be useful in analysis. Fung [3], Chapter 6, 
lists many such possibilities by indicating the appropriate form of 
the loading function, and his bibliography cites much of the literature 
in this area. Both isotropic and anisotropic forms are given, the 
latter deriving from Koiter's work [7], and various examples are cited 
that provide for kinematic and isotropic hardening, and for a 
Bauschinger effect. Indeed, Fung provides an informative overview of 
this aspect of plasticity without becoming involved in the details 
of individual relations or solution methods. 

It is useful here, however, to examine such details for a few 
situations. In so doing, we are altogether mindful of the dual needs 
for computational convenience and physical sensibility. Thus, for 
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example, we do not discuss loading functions that reflect a 
Bauschinger effect; the experimentation required to infer the form 
of such a function alone is prohibitive in the absence of compelling 
reasons for such work. Still, the procedures outlined are easily 
adapted to such representations when warranted. 

Anisotropic Relations: Following Koiter's generalization [7], 

we consider not one but a set of n load functions 

f (a. .,ef?]w^ p b - 0 a = l,2,...,n (22) 

The condition f =0 describes a segmented surface in stress space 
comprising n mutually exclusive segments. These pieces of course 
are contiguous and, since the stress vector corresponding to any 
point in the physical domain must be single-valued, yield proceeds by 
activating only one load function at any instant of time. Thus, if 


then 



0 


f < 0 a = 1,2,..., 3-1, 3+1, ...,n 

Of course, different f may be operative at different points of the 
physical domain at a given time, or at the same point but at different 
times. If for convenience we let 


f (o..,e^, IV^) = <f> (a..)-i|> (e. (p) ,W tp) ) 
ij ij 1 Y a v \y r a lj 1 


( 23 ) 
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it is a direct calculation to arrive at the flow rule 


11 


n 


= E 

ct=l 


2p( P) (34 /8o )o 

a a mn mn 


6 v 


( 24 ) 


as an analogue to (11) . In (24) ■ is a set of plastic moduli 

analogous to in (11). Note that since flow is proceeding on 

only the 0th surface, the summation in (24) is not a collection of 
processes that occur simultaneously, but of those only one of which 
may occur at any time and position in the physical domain; (24), 
however, represents the totality of processes that may go forward 
at all times and points. 

The simplification (23) is non-essential and an analogue to 
(8b) may be written. It would involve replacement of f by f and a 
summation, as in (24), and need not be written here explicitly. To 
either such a result or (24) may be added the necessarily linear 
relations between the elastic strain increments and the stress increments; 
the result assumes the form (15) which is necessary for the remaining 
development of the theory. It should be noted that the elastic relations 
may be fully anisotropic with no conceptual difficulty. There is 
certainly the potential tedium of inverting such equations to arrive at 
the analogue to (16), but, to the extent that we envisage extensive 
use of the computer in solving problems, no practical impediment to an 
anisotropic theory of elasto-plastic flow is apparent. 
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Isotropy and Incompressibility : Up to this point, we have allowed 

the function (j>( 0 ^) (or the functions ^(a^) to be arbitrary. 
Experientially, however, it is clear that specific functions are very 
attractive from the standpoints of computational convenience and 
physical sensibility. In particular, we may now require flow to 
be both isotropic and incompressible by restricting <j> to depend 
only on the stress deviator 


S. . = a. . 

il il 




and, more particularly, its invariants 


J 2 = S . . S . . /2 


il li 
ij S jk~ki' 


J 3 = S^.S.^./S 


so that, for example 

<j> = <J> (J 2 > J 3 ) 

This form is appropriate to the observations made by Osgood [8] and 
discussed by Drucker [9] wherein 

* = (j| - C9/4)J3> 1/6 


(25) 


(25a) 


is suggested as an equivalent stress. The resulting equations pertinent 
to (25) were given in some detail in [6] and need not be repeated here. 
It is a straightforward matter to use the specific form (25a) in these 
formulae, althouth we have not pursued such an effort. 

Alternatively, one may revert to the simpler and more popular form 
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associated with the name of von Mises. In this case 


' <f> = /(2J 2 /3) (26) 

and the equivalent stress x g ^ is indistinguishable from the octahedral 
shear stress x . Using an octahedral plastic strain increment in 
place of the equivalent quantity defined in (9) , we find that if 


then 


4x /Sc = 2 V tp) 

O 0 0 


p ^ = 3 V tP) 

o eq 


It is then easy to show [5] that the flow rule is 


2vSc ij ‘ Sa if^^ Sa a 6 ij * ^0 P) >Vk£ 5, V /3T o 


and its inverse is 

,/2u = - S ij S k£ «e fce /[3T*(l*P 0 tP) /p) ] 

The extended form of Navier's equations is, in the absence of body forces. 


[l/Cl-2v)]6u._. it 6u iJj -2fS..S k£ 6u i y[ 3 ^tlxp o tl,) /»)) i 


= 0 


(27) 


and it is again useful to note the quasilinear form of this diffemetial 
equation: The dependent variable 5u^ appears linearly in (27), while 

the nonlinearities of the problem are confined to the coefficients of 
the last term. The operational result of this important feature is 
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discussed below as well as in separate publications dealing with special 
aspects [5,10]. Finally, we note that, as y ^/v 0 (perfect 

plasticity), the flow rate in the form shown loses meaning; however, 
both the inverse and the differential equation (27) go smoothly to the 


limit and exist and are tractable. 
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FINITE DEFORMATIONS 


of the deformations and the implied assumption a../y^ <<: 1 

ij eq 


Up to this point we have developed the theory of elasto-plastic 
flow wholly within the context of infinitesimal strains, that is, for 
e^. sufficiently small for (17) to apply. As Rice [11] and Hill [12] 
have noted, however, there may be instances in which the magnitude 

are 

such that this framework is inappropriate. For these and other 
reasons, the theory was extended to the case of finite deformations. 
The development is presented in considerable detail elsewhere [13], 
and we give below just its main points.* 

The total deformation is described in a fixed reference frame 
by the mapping 

x 1 = x 1 (X J ,t) , | x 1 , j | i 0 


(28) 


In (28), x 1 are spatial coordinates of material particles comprising 

domain B at time t -> t which were located in B at t . We 

o oo 

differentiate (28) with respect to time to obtain a velocity field 
v 1 = v 1 (x' , ,t) and compute the velocity gradient 


v. . = d. . + u). . 
nj ij ij 


*For this purpose, we use in this section general tensor notation, 
covariant (contravariant) character is denoted by subscript (super- 
script) indices; a comma (semicolon) indicates partial (covariant) 
differentiation; repeated indices in subscript-superscript pairs 
implies summation over the range 1,2,3; x 1 and X* are coordinates in a 
single, fixed, orthogonal curvilinear system yi; g. . is the metric 
tensor of y*; and <$i is the Kronecker delta. 1 - ) 



27 


where d. . is the (symmetric") deformation rate tensor and u). . is the 
ij ij 

(skew-symmetric) spin tensor. These tensors are of particular 
significance in that they both provide a full description of motion, 
and are linearly dependent upon v 1 . 

We next invoke the assumptions of the existence of a load function 
and Drucker's hypothesis. In his original work [13], Osias used a 
loading function of the form of (25) — in the context of (3) — but it 
seems clear that far more general forms are altogether acceptable. 

The statement of Drucker's hypothesis is written as* 

o^d^P > 0 (29) 

ij 

i i 

where a J is the Jaumann stress rate given by 


-ij *ij i kj j ik 
a J = a J +a w 


= 3a^^/3t+(9a^/3x^)v^ 

and d^j is the plastic component of the deformation rate. It may 
be shown that, independent of any rotation, (29) requires that work 
done by the external agency in producing d:^ is nonnegative, con- 
sistent with the sense of the original hypothesis. The Jaumann rate 
is used in (29) to provide an objective measure of the change in stress 
viewed from a frame rotating with the material, as follows from the 
work of Prager [14]. 


*Note that (e) , (p) , and eq are not indices. 
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With these bases, Osias derives a flow rule of the form (15); 
the details, however, are such that explicit notation is in order. 

The flow rule is 

- B i jk £ S “ (30) 

B ijk£ = ^ ik g j£ + g i£ g jk )/2-[W(Uv)]g i;j g k£ 

+ (3<|)/3c k ^) 

Y = (p/^?)U/[(94>/3cT ij )a ij )} 
and the inverse is 

a 1J /2y = D l:iM d u (31) 

D ijk£ = Cg^VVV^/tl^vJJg^g" 

1/y + rS6/8c m j (5)/Sc rs jg mr g riS 

and g 1 ** is the associated metric. Note that the elastic deformation 
rates and the Jaumann stress rate are presumed to be related in precisely 
the same fashion as are infinitesimal strains and stresses — or their 
increments — in classical isotropic elasticity. To that extent, 

(30) and (31) are tailored to materials for which a shear modulus and 
and a Poisson's ratio may be defined in the elastic range. It should 
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also be noted that (30) and (31) do not involve strains per se; such 
quantities are computed a posteriori in analysis. 

In requiring equilibrium, Osias adopts the rate or increment 
viewpoint which has proven so useful. He observes that the total 
force acting on a body must vanish, i.e.. 



0 


where t 1 are surface tractions, and requires that the time rate of 
change of this integral also vanishes. Hence 



c,i: ' v \kV dS 


0 


so that the traction rate is 


•i .*ii ki j . 
t = (a J -a v J >k )rij 


(32) 


(33) 


and the equations of equilibrium are 






. -a 


ki 


;j ;k 


= u 


(34) 


with no body forces acting, of course. In the interior of B, (34) 
apply. On the boundary S, however, there is some choice. One may 
specify velocities in the usual way; one may specify traction rates 
(33) when, for example, a pressurization is known; or one may specify 
a force in terms of an integral as on the left of (32) when such a 
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quantity is known over some portion of the boundary S. 

Combining (34) and (31), along with the definition of deformation 
rates in terms of velocity gradients gives the final equations for the 
latter as 




it kj 


[ D v j;d ;i * 8 o 'v,., 


i; j l 


(35) 


,, 3k .3 £k. .k.3 1 1 A.£ 13, r mn _ 

+ (1/2)[6„((S g )-<5 <5 J g -6 6 g J ][o v. »] . = 0 

v J L n v m 6 m 6 J m n 6 n m 6 JL 3 ;£ j ;i 


Inspection of (31) and (35) shows the latter to be quasilinear in 
the velocity v^ so that, aside from increased complexity, (35) is 
operationally of the same form as (20) or (27) . 

As a result we have a form of Navier's equations further 
extended to the case of elasto-plastic flow of bodies permitting 
finite deformations, together with boundary conditions appropriate to 
their integration. The various comments in preceding sections of 
this report concerning the utility of this sort of formulation 
pertain equally here. While we have not examined such matters as 
integral theorems, we have used a planar counterpart of (35) to solve 
problems, and these results are described fully elsewhere [13]. 
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REDUCTIONS TO 2 -SPACE 


While derivation of the theory is easily prosecuted using indicial 
notation, it is useful to reduce the final relations to a form akin to 
that employed when actually solving problems. Accordingly, we devote 
this section to the governing equations for a series of situations in 
2-space, that is, in more familiar extended representations. We 
perform this reduction only for Mises or J 2 theory where <p is given 
explicitly by (26). Forms appropriate to (25) are given in [6]. 


Torsion: There are two cases of torsion, one for prismatic 

bars or bars of constant cross-sectional shape, and the other for 
axisymmetric rods whose shape is otherwise variable. The first may 
be described in rectangular cartesian coordinates (x,y,z), and the 
dependent variables are the corresponding displacement increments 
(Su,<5v,5w). It may be shown that, if the z-axis lies along the 
center of rotation, 


<5u = 

-yz<Stf 

6v = 

xz<5t? 

6w = 

6w(x,y) 


where 6t? is an increment in the twist per unit bar length. It follows 
that the stresses are given by 

0 =0=0=0 =0 



y z 
= T x (x,y) 

= T y (x,y) 


xy 
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and that equilibrium requires 

9<St /3x + 3<St /3y = 0 

x y 

We define 

e - 

and observe that 



to write the flow rule 

pSw, x = yy<5t? + Cl+eT 2 x /J 2 ^r x + (BT x T y /J 2 )uT y 

y6w, y = -xySi? + ( 6T y T x / J 2) 6t x + (l + 3x y 2 /J 2 ) T y 
and its inverse 


[(l+B)/y]Sx x = [(1+Bx 2 /J 2 ) (6w, x -y6t>)-(BT x T y /J 2 )(5w,^+x6^)] 


[(l + B)/y]6x y = [-(6T y T x /J 2 )(6w, x -y6t?) + (l + eT2/J 2 )(6w, y+ x6^)] 

If we apply equilibrium to the inverse of the flow rule, we obtain 


(1+Bx 2 /J 2 ) 6 w, -2(Bx x /J 2 )6w, + (1+Bx 2 /J 2 ) Sw , 

y xx x. y xy x yy 


+ a^5w, + a 2 6w, = (yai-xa 2 )6t? 

x y 


where 


«1 = a + 3){ [(l+8t 2 /J 2 )/(l+g)], x - [CBT x T y /J 2 )/(l+3)], y } 
a 2 = (l + B){-[(BT y T x /J 2 )/Cl+B)], x + [(1+Bx 2 /J 2 )/(1 + B)], y } 


(36) 


( 37 ) 
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Boundary conditions for (37) are 

6t cosit + 6 t sinii) = 0 

x y 

where \p is the angle between the unit outward normal and the x-axis, 
and the stress increments relate to 6w and its derivatives through 
(36). Note that, where no flow occurs, 8 = 0 and (37) reduces to the 
familiar 

V 2 8w = 0 

The alternative formulation derives from elimination of 6w from 
the flow rule, having introduced an extended Prandtl stress function as 


t x = 9?/9y 


t y = -9c/9x 


with C = c(x,y). The differential equation is 


(1 + Bd x /J2)^, xx + 2(B?, x C, y /J 2 )6C, xy + (l+8C, y /J 2 )6c, yy 


+ a 3 6c, + a 4<5 C, v = -2y<5t? 

a y 


(38) 


where 


= VH 


+ 


a 3 = (l+BC, x /J 2 ) , x + /J2)» y 

a 4 = (65, y C, x /J 2 ), x + (1 + 8C » y /d 2 ) > y 


and 6^ = 0 on the boundary of the domain so long as it singly connected. 
This last restriction renders (38) less useful than (37) , for multiply 
connected domains are excluded from study. Still, the quasilinear form 
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of (38) is evident, as is true for (37), and solution procedures may 
be developed to exploit this feature. 

Torsion of axisymmetric bodies is described in a similar fashion. 
In a system of cylindrical coordinates, the only active displacement 
is in the circumferential direction and we denote it by v(r,z) . The 
only active stresses are 


a . = x 
r0 r 


a = t 
6z z 


and 


T 2 2 
Jo = T + T 

r 7 


Using the same definition of 8, the flow rule becomes 


y6, z = (1 +Bt z /J 2 )<5t z + (Bx^/Jz) 6x r 

v(<5v, r -Sv/r) = (6T r x z /J 2 )'5T z + (1+3 t£/J 2 ) 6x r 


and its inverse is 

[(l+0)/y]6x z = (H-J3 t*/J 2 )5v, z - (Sx^/Jz) (6v, r >6v/r) 

[(l+3)/y]6x r -(3x r x z /J 2 )6v, z - (l+gx^/J 2 ) (6v, r >6v/r) 


(39) 


If we apply equilibrium 

96x /9r + 2Sx /r + 96x /9z = 0 
r r z 

we have 

(l+6x z /J 2 )6v, rr -2(8x r x z /J 2 )6v, rz + (1+0 t^/J 2 ) 6v, zz 
+ a'(6v, r -6v/r) + a 2 <Sv, z = 


0 


( 40 ) 
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where 

< = (1+6) { [(l+6T^/J 2 )/(l+6)], r - [C6T z T r /J 2 )/(l+B)], z }+(l+6T^/J 2 )/r 
« 2 = a+6){-[(SVz /J2)/Cl+8)3j r + [C 1+ ^J/J2)/(l+6)], z }-(eT r T z /J 2 )/r 

Boundary conditions for (40) reduce to 

(St^cosiJ) + 6x z sin^ = 0 

where $ is the angle between the unit outward normal and the r-axis 
(i.e., the plane z=0) , and the stress increments relate to 6v and its 
derivatives through (39). Note that, for purely elastic behavior, (40) 
reduces to 

6v, rr +(1/r3<5v, r " + 6v » zz = 0 

which of course is the standard elastic result.* 

Alternatively we may introduce the stress function t(r,z) through 

r 2 r = -C , , r 2 r = S , 

r ’ z z r 

and eliminate reference to <5v in the flow rule. The result is easily 
written terms of the ratios 

C r = ?, r /|V5| , ? z = G, Z /|VC| 

with |V ? | = /(£, 2 + 5 ,^) 

*A slightly altered form of both this relation and (40) is obtained if 
the dependent variable <5v is replaced by an angle of twist 6t?(r,z) = 
6v(r,z)/r. 
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as 


(l+BC;)[6c >rr -(3/r)6s, r ] + 


2 Vz 5 Wz +( 1 ^z )6? >zz 


with 


+ + 0 ‘ 4 <S ?> = 0 

a 3 = C 1+ ^’r + 

<*4 = (B? z ? r ), r + (1 + BC Z ), Z - (3/r)B? z C r 


(41) 


The boundary condition on (41) is of course constancy of 6£, as in 
elasticity. In both (40) and (41) the quasilinear nature of the 
differential equation is evident as noted for other cases in this 
section. 


Axisymmetric Extension: Still operating in cylindrical coordinates 

(r,0,z), the non-zero displacements are observed to be 

u = u(r,z) 
w = w(r , z) 

and the strains are 

e = 9u/9r 
r 

£ e = u/r 

e = 9w/9z 
z 

Y rz = 9u/9z + 9w/9r 

to which there are corresponding stresses a , a., a , and x . The 

To Z T Z 

flow rule is conveniently written in terms of normalized deviators 



37 - 


= (2a 


r“V<y/ ( 3 / 3 t q ) 

( 2 a e~ a z ~ a T )/( 3 / 3 r ) 

cr 

(2CT z" a r -^ 0 )/(3/3T ) 


rz 


where 


We £ hen have 




T 2 
0 


C 2 / 9 )(a 2 + a 2 + a 2 -a „ 

r 6 2 6 z~ a z a r' cr - cr ~ + 3T 2 


r e 


3t ) 


rz^ 


2]j<$ 


e r = fVa + v) +6s 2 J(Sa + . 

r r L -V /(l+v)+g S s To r 

r qJ 5 e + f-v/(i +v ) + g s ^ s ^ ^ 


+ 2s s 

r rz r z 


2jJ<Se 0 = f'V/(i +v) + 


8S S S r ]5 Vfl/(l + v) +6s 2 J(Sa 

0 0 f v/ ( 1 + ^ + Bs 6 s]Sa 

0 2 2 


+ 2 s s 6 t 
9 rz rz 


2y6 


2 * ( - v/ n^8s z s r]6vr . v • 

6S Z S 9J 5 Vtl/aev)e es 2 JSc| 


+ 2S 7 S rr^ T 

z rz rz 


Z J Z 


W5y rz = 6s T^ s „<Sa + 
aS the flow rule. 'The 


rzVV g s s n <$a_ + f 


rz e Oa 0 + es r2 S z 6a 2+Cl+ 3 s 2^ 6T 

averse is written 


rz 


a " d the La, " e ' instant A. „ here “ Slng V °“ nS ' S " odul us E 

A = 2vy/( i _ 2x; - ) 

7 » £ = 


2(J+v)j 
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and we have 

Sa r /E = [ (X+2y) /E - s*/(l+l/8) ]6e r + [X/E-s r s 0 / (1+1/8) ]6e 0 

+ [X/E-s r s z / (1+1/8) ]6e z - [s r s rz / (1+1/8) ]6y rz 
6a Q /E = [X/E-s 0 s r /(l+l/8)]6e r +[(X+2y)/E-s 0 /(l+l/8)]5£ 6 

+ [X/E-s e s z /(l+l/8)]6e z -[s 0 s rZ /(l+l/8)]6Y rz 

(42j 

6a z /E = [X/E-s z s r /(l+l/8)]6e r +[X/E-s z S 0 /(l+l/8)]6e 0 

+ [(X+2 U )/E-s z /(l+l/8)]6e z -[s z S rz /(l+l/8)]6Y rz 
<St rz /E = -[s rz s r /(l+l/6)]6e r -[s rz s 0 /(l+l/8)]6e 0 

- [s rz S z /(l+l/8)]<Se z +[y/E- S 2 z /(l+l/8)]6Y re 

Navier's equations are obtained by combining (42) with the equations 
of equilibrium, viz. 

86a /8r + (6a -6a Q )/r + 86 t / 8z = 0 

T I* a X* Z 

86t /8r + 6 t /r + 86a /8z = 0 

rz rz z 

having expressed the strain increments in terms of displacement increments 
in the usual manner. Boundary conditions are given in terms of 6u, 6w, or 
tractions which are linear combinations of (42) . The result is more 
involved than (40) or (41) , say, but the operational nature is very 


much the same. 
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Planar Behavior: It is useful to consider two types of problem 

under this heading, plane strain, and combined plane stress and 
(Kirchhoff) bending. The first may be established by formally requiring 


e z = 3w/3z = 0 

Y yz = 3v/3z + 3w/3y = 0 

Y xz = 3u/3z + 3w/3x = 0 

where the domain of interest lies in the x-y plane. We then operate 
on (27) and its associated relations to develop the requisite equations. 
From (17) we have 


<5 e = 3<5u/3x 

x 

<$ e y = 36v/3y 

6y x y = 3<5u/3y + 3<5v/3x 

and the flow rule is again written in terms of normalized deviators: 

s = (2 a -a -a )/(3/3x ) 

x x y z o 

s = (2 a -o -o )/(3/3x ) 

y y z x' o 

s = (2o -a -o )/(3/3x ) 

z z x y J o 

s = x /(/3x ) 

xy xy' o J 

with 


x 


2 

0 


2 2 2 

(2/9) (o +o +a -oo-oo-oo +3x 
v x y z yz zx xy 



We note that the transverse stress o z accumulates from observing the 


0 = 


(l+gs z )6a z +(-v+Bs z s x )6a x 

+ (-v + es z s y )6a y+ (2gs zSxy )6x xy 


constraint 
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which, together with <St = 6t 

& xz yz 

The flow rule then becomes 


0, insures the plane strain condition. 


2y Cl + Bs^)6e x 


[1-v+6(s|+2vs x s z +s 2 z )]6o x 


+[-v+g(s s -2vs^)16a 
L x y z J 1 y 

+2[g(s +vs )s ]6x 
L v x z J xy J xy 


2}j(1+Ss )6e 
z y 


[-v+g(s y S x -2vs z )]Sc x 

2 2 

+[l-v+g(s +2vs s +s )]6a 
y y z z J y 


+2[g(s +vs )s ]6 t 
L y z' xy xy 


p(1+es z )6Y xy 


res (s +vs )]<5cr 
L xy x z J 1 x 


+ [3s (s +vs )l6a 
1 xy y z J y 

+{l+2g[s 2 + (1+v) s 2 ] }<5x 

xy z J xy 


with the inverse 

6a x /E = [CA + 2y)/E-s x /Cl+l/g)]6e x +[X/E-s x s y /Cl + l/g)]6c y 
-[ S x S xy /(1+1/3)]6Y xy 

Sa y /E = [A/E-S y S x /Cl + l/g)]6E x+ [(A + 2y)/E-s 2 /(l + l/e)]6 Ey 


-[s s / (1+A/g) ]6y 
1 y xy J 1 xy 

6x xy /E =-[s xy s x /(l + l/6)]6 Ex -[s xy S y /(l + l/g)]6e y 
+ [y/E-s^y/ Cl+1/6) 3 « Txy 
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which are then combined with the equilibrium equations 

36a /3x + 36 t /3y = 0 

x' xy 3 

36t /3x + 36a /3y = 0 

xy y 

to produce the final result . 

The combined bending-stretching equations are derived by using 
the analogue to to the theorem of minimum potential energy, outlined 
above. We consider the displacements and their increments to be 
functions of the coordinates (x,y) only, but write the strain 
increments as 

6e^ = 36u/3x - z3 2 6w/3x 2 

6e^ = 36v/3y - z3 2 6w/3y 2 (43) 

6y^ = 36u/3y + 36v/3x - 2z3 2 w/3x3y 

and insert these' into the functional II under the condition that 

t =-r=a = 0. Jones [4] has worked this out in full detail, 

xz yz z L J 

so we shall not repeat his results here. It may be noted that he 
recovers a coupled theory, and this coupling occurs both in the 
field equations and in the boundary conditions. The former are in 
terms of the three displacement increments (6u,6v,6w) and are of 
course quasilinear. The latter are in terms either of these three 
quantitites or resultant increments of force (for in-plane loading) 
and of moment and equivalent shear (for out-of-plane loading) . 
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He also demonstrates appropriate 'corner conditions', in terms of 
twisting moment and comer angle. Through solutions to problems, he 
shows that, although the strains remain linearly distributed through 
the plate's thickness per (43), the stresses become nonlinear as 
yield progresses. He also makes clear the manner in which bending 
and stretching couple in the presence of plastic flow, although. the 
two events are distinct for purely elastic behavior. We do not 
reproduce his findings here both because they are accessible elsewhere 
and because their form is evident from the foregoing developments. 
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SOME ELEMENTARY SOLUTIONS 

It may be useful, in order to see more clearly the workings of 
the theory, to consider some elementary examples. We therefore look 
in this section at torsion of an axisymmetric rod of constant radius 
R, and at bending and extension of a rectangular bar. For the first 
case, the circumferential displacement increment is simply 

Sv = rz6t? 

where, as in the last section, 0 is the twist per unit rod length. 
The flow rule is simply 

prSi? = (1+6) St z 

having used the notation of the last section. Since the only active 
stress is t and 

6 = U /„ 0 (P) 

where y Q = y ^ (x ) and t q = /(2/3)x z , the flow rate may be easily 
integrated. The result is 

yri? = t z + 2ye ( ^ P ' ) [/(2/3 )t z ] 

and is a given function of or, as implied, /(2/3 )t z . Next, 

we observe that the twisting moment is simply 



(44) 


o 
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so that moment and twist may be related. Let. us suppose, for example, 
that 

' (rt - i 

J 1/n 

T > T, 


0 


T < 

o Z 


e * (T o /t £ _1) 


where e*, t^, and n < 1 are constants. Then if 


T = (2/3}t £ 


we have 


y#r = 


t + 2y£*(r /t -1) 
z v z p ' 


1/n 


0 < r < x /yt? 
P 


x /y < r < R 
P 


Thus the elastic relation is recovered for Rt? * t /y and 

P 


(45) 


= yJt? , 


J = ttR 4 /2 


but as yield progresses, (45) obtains. In that event (45) must be 

evaluated for as a function of r to find M. . As an illustration, 
z t - 

let n = 1; we have then 


„ T 0 + 8e*/3R 

M. = yJ , ~ — 

t 1 + 2yc*/x 


which gives a moment-twist relationship that is asymptotic to the 
implied bilinear 'curve'. The case of perfect plasticity may be 
represented by e* -* °°, in which case 

M , . = (2tt/3)t R 3 

t-lim P 
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is a limiting moment that may be applied to the rod. 

For rectangular beams, Jones [4] derives the governing equations 

with the aid of the analogue to the theorem of minimum potential 

energy. The beam occupies the region 0 < x < L, -y /2 < y < y Q /2, 

-z /2 < z < z / 2: the end at x = 0 is cantilevered; the end at x = L 
o o 

is subject to increments in an axial extension force 6 F, a shear 
force 6 Q, .and a bending moment 6 M; the length of the beam is loaded 
transversely by Sq. Since the flow rule is 

6 a = (1/E + l/3p^)6a 

where subscripts on 6 e and 6 a are suppressed, and the inverse 
(including transverse straining) is 

6 a = E [1- l/(l+3y.^ P - > /E)]6e 

he takes the strain-displacement equations in the form 

6 e = d 6 u/dx-yd 2 6 v/dx 2 

and computes the functional 31. Operating to find a minimum value, 
he obtains the equations 

d (Cid 6 u/dx)/dx - d (C 2 d 2 6 v/dx 2 )/dx = 0 

( 46 ) 

d 2 (C 3 d 2 6 v/dx 2 ) /dx 2 - d 2 (C 2 d 6 u/dx)dx 2 = 6 q 



- 46 


where 


Ci = z 


I 

-V 2 

r 


C 2 = z Q | yE "dy 

V 2 


r 


C 3 = z o | y 2 E'dy 


and = E where no yield occurs, otherwise 


E " = E[1-1/(1 + 3^ P) /E)] 

The boundary conditions (at x = L) are 

Cid 6 u/dx - C 2 d 2 5v/dx 2 = 6 F or 6 u 

C 3 d 2 6 v/dx 2 - C 2 d 6 u/dx = 6 M or dSv/dx 
d(C 3 d 2 < 5 v/dx 2 )/dx - d(C 2 d 6 u/dx) dx = - 6 Q or 6 v 


specified 
specified (47) 
specified 


It is seen that C 2 couples bending and stretching in both the 
differential equations and the boundary conditions. Furthermore, 
all of the coefficients (Ci, C 2 , C 3 ) are problem dependent so that 
solution to a given problem may not procede in a simple manner. 

If the bar is loaded in simple tension, then it is easy to show 
that the original stress-strain relationship is recovered, albeit in 



47 


uniaxial quantitites. On the other hand, if the bar is put into pure 
bending, matters are more complex in that there is an elastic core 
sandwiched between two elasto-plastic outer layers. The strain 
distribution remains linear through the beam's thickness, i.e., with 
respect to y, but the bending stress distribution becomes nonlinear 
as a reflection of the stress -strain curve. 

For the particular case of pure bending, is of no interest; 

C2 is null owing to the symmetry in E'; C3 has different values in 
the core and outer layers, the latter depending on the shape of 
the curve; the increment in curvature is given by 

6k = d 2 6v/dx 2 

and the (incremental) moment -curvature relation is 

C 3 6k = 6M 

which satisfies both ( 46 ) and ( 47 ). So long as a is everywhere 
less than the yield stress Y, we have 

(Ey 3 /z /12)6k = 6M 

and the deflection increment is 

6v = 6 [6M/(Ey^Z Q ) ]x 2 

which are the familiar elastic results, and they obtain for 

M < My = Yo Z o Y / 6 

As the moment M exceeds My, however, we need to account for the stress- 
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strain relationship. Suppose that 


.(P) _ 


e u ' = e* (t /t 0 -1) 
o o l 


as before, except that the exponent is by implication set to unity. 
Then 

6a = E { 1 - l/[l+Y/(/2Ee)]}6e 
and, since 6e = ySx, we have 


6M = 2z 


y 2 v 2 

0 J Ey 2 6Kdy + 2z q j E{1-1/ [1+Y/ (/2Ee*) ] }y 2 6Kdy 


V 2 


- CEy 2 z Q /12) ( fy p /y 0 ) 3 * U- 1/ [ 1 *y/(,'2Ee*) ] } [l-(y p /y 0 ) 3 ]j 6 k 

where y^ = 2 Y/Ek. Observing that y^ depends upon k, we integrate this 
last expression to obtain the final relationship. Setting 

k* = Ey Q x/2Y 

8* = l/[l+Y/(/2Ee*)] 

we have 


M/My = - g*/2fc* 2 + 38*/2 
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the additive constant of integration being adjusted to connect this 
result with that for purely elastic behavior. At k* = 1, M/My = 1; as 
k* becomes large M/My -> (1-£*)k* + 3g*/2. In addition, perfect 
- plasticity is characterized by 8* -*■ 1 so that we observe the limiting 
situation of M/My -*■ 3/2. For perfect, plasticity, it is easy to show 
that 



so long as 1 < M/My < 3/2. 

Thus we have solution to some elementary problems. Slightly more 
complicated problems may of course be devised, but the algebra becomes 
far more difficult and it is not in our interest here to pursue the 
matter. The point remains, however, that a variety of approaches is 
available to solve problems and that, while we may develop explicit 
formulae for simple cases, we are always in a position to observe the 
process of response from purely elastic to combined elasto-plastic 
flow. We may also allow excitation to reverse and examine the process 
whereby residual stresses develop. 
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SOLVING MORE COMPLEX PROBLEMS 

To this point in the development, it should be clear that the 
theory of elasto-plastic flow is fully described by a set of one or 
more differential equations, and boundary and initial conditions. 

The governing differential equations are typically* written in terms 
of displacement increments; they are quasilinear with variable 
coefficients which depend upon the current stress state and given 
material properties and thereby reflect the nonlinear aspect of the 
process. The boundary conditions are in terms of displacement increments 
or traction increments which depend linearly on the gradients of 
displacement increments. Initial conditions perforce are to be given 
in terms of stress and strain fields. 

Problem solving will usually entail numerical procedures and use 
of a large, high-speed digital computer. The issue then becomes , 

identifying an appropriate method. If we look at the kind of problem 
to be treated, as summarized above, we observe that the equations for 
incremental quantities are similar to those for ordinary elasticity, 
but with two important exceptions. The first is that the apparent 
material properties, as described by the coefficients, are inhomogeneous 
and effectively anisotropic. The second is that distinction must be 

*Altemate formulation in terms of stress function(s) is permissible 
but, on the whole, has not proven to be particularly useful. Other 
than the one example shown above, we do not pursue this form further. 



51 


made between loading and therefore continuing elasto-plastic flow, and 
unloading and thus only elastic behavior. 

These exceptions give no especial impediment to problem solving 
when a computer is used, for these devices are easily programmed to 
handle what are, after all, merely complex details. Accordingly, the 
question of selecting a method appropriate to elasto-plastic flow is 
no more difficult than settling on a numerical technique for ordinary 
elasticity. 

Accordingly, we may choose from a more or less standardized list: 
finite differences, finite elements, and so on. In the specific problems 
we have solved and released for publication [15-19], we have preferred 
finite elements of the simplest kind, largely for reasons of convenience 
and economy. These considerations have proven sensible in dealing with 
more complex situations of bending [4] and finite deformations [13]. 

We have also formulated an approach using the boundary-integral equation 
method [5] , and Riccardella has implemented this technique for the 
case of plane strain [20] . The overall algorithm has been described 
as well, in a context that transcends specific method [21,22,10]. In 
most methods of practical interest, differential equations of the 
form (27), but reduced to a specific situation, are then reduced to 
algebraic equations which we represent as 

= { 6 1 } 


[A] (6u) 


(48) 
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with boundary conditions being specified in terms of certain elements 
of (Su) and/or {St}. In (48), [A] is a matrix of coefficients that 
replace the differential operator in (27) , {Su} is a vector of 
displacement increments at points distributed throughout the domain 
of interest, and {St} is a vector of traction or force increments 
at these points such that, except for specified boundary locations 
{St} is largely null. Note that (48) obtains for both finite 
differences and finite elements, and that a parallel form is sufficient 
to represent the equations for the boundary integral technique. 

Program Outline *: Typically, a computer program for elasto- 

plastic flow in two spatial variables comprises an executive routine 

followed by one that sets up data and initializes various arrays, and 

then by several operational routines . The initializing routine reads 

in geometric information, loading data, and material properties. The 

load data may consist of tractions and/or displacements, in either 

of two orthogonal directions. The material data include two elastic 

constants and the relation between and x . This is read in as a 

o o 

table of coordinates of a sequence of points along the octahedral 
stress - plastic octahedral strain curve; spacing between the points 


*This discussion is taken largely from [10] ; we appreciate the 
cooperation of that publication's editor in releasing the text for 
presentation here . 
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is roughly proportional to the curve's radius of curvature. It is also 
useful to provide an artificial extension to the curve, and we assume 
indefinite continuation of the slope on logarithmic scales, as evaluated 
near the end of the input data. Alternatively, one could employ a 
simple straight line extension. Use of this technique precludes 
inadvertant loss of a run because the stress-strain curve was out-run; 
the analyst is enabled to decide where to cut off his results after 
they are in hand. 

The first load step is necessarily elastic, and no difficulty 
is encountered in performing it for an arbitrary level of loading. 

The matrix [A] is assembled, the vector {6t> is computed together with 
whatever modifications in [A] are required due to imposition of boundary 
conditions, and the resulting equations are solved for {<5u}. Deformation 
and stress increments are determined. Excitation may then be scaled so 
that the most highly stressed point in the domain just exceeds the 
proportional limit, e.g., T c / T £' = 1.00010. 

Before beginning the next load step, certain operations are 
necessary, so that the initiation of yield may be taken into account. 

The table of points describing the stress-strain curve is scanned, 
beginning with that at the proportional limit (i.e., (0, in the 
table described above) . The scan continues until a point is found 
whose coordinates exceed that of the point in the material; the next 
highest point is then selected for further computation. The slope of 
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the curve (2y^^) is then taken to be that of a straight line connecting 
this point on the curve and the material point. 

This procedure accomplishes several things. It is clear that the 
value of thus determined is a reasonable approximation to an average 

over the increment of loading to come. Proximity of the material 
point to the point in the curve will not exert undue influence on the 
value of which, incidentally, is necessarily positive. If the 

material point happens to be off the curve, it will be steered back 
on; both overshoot and undershoot are corrected. Furthermore, 
extensive experience has demonstrated that this procedure is inherently 
self- corre ct ing . 

At the position where yield has just begun, no significant plastic 
strain or strain energy has occurred, and no values of these quantities 
are stored. The calculation proceeds with whatever further data 
reduction and output preparation are appropriate. In our finite 
element programs, for example, we find strains, stresses, and several 
of their invariants. As is evident from (37), for example, at least 
the stresses are needed to prepare for the next load increment. 

The second load step follows the same pattern used for the first, 
with two exceptions. Assembly of [A] is modified to account for yield 
through introduction of suitable material parameters where yield has 
begun: Also, after finding {Su}, the same considerations obtain in 

determining the corresponding stress increments. 
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We have preferred to reassemble [A] fully for each load step, 
rather than to modify a retained and originally elastic version of 
this array. While there is no unanimity over this point among 
practitioners in this sort of analysis, it is pertinent to note certain 
bases of our preference. Fresh reassembly of [A] for each load step 
is not, in our experience, costly in time; the core storage allocated 
to [A] may be overwritten for temporary purposes once {6u} is determined; 
the more complex bookkeeping needed to modify [A] from step to step 
is avoided and the overall code is thereby more compact; and this choice 
directly admits multi-cycle loading in that different types of boundary 
conditions may be imposed from one step to the next. 

Subsequent steps follow the pattern outlined above, and the 
computation nominally continues to whatever termination point has been 
set. Because, however, unloading either locally or globally may occur 
during any load step, we insert an additional procedure prior to finaliza- 
tion of the results for each step. Once increments in displacements, 
strains, and stresses are found, the solution is regarded as having 
only a candidate status. The program checks to determine whether 
the candidate solution implies an increase in t q at all points currently 
yielding. If so, the candidate solution is accepted and the program 
moves on to completion of the load step. 

If not, then [A] is reassembled with elastic constants in those 
positions corresponding to material points where elastic unloading 
seems to be in process. The degree of yield at which unloading is 
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(tentatively) first detected is stored, and an iteration within the 
increment is begun until there is internal consistency between the 
presumed and actual behavior at each material point. Once such 
consistency is achieved, whether or not iteration is required, the 
increment in plastic octahedral strain is computed via 

Se (p) . 2 m Cp) 6t 

O 0 0 

and the plastic work is obtained through simple trapezoidal quadrature 
at those points where yield continues. 

Where unloading has occurred, t q is traced so that when it 
again reaches the level where unloading began, the elasto-plastic flow 
rule is again used for computation of [A], and plastic strain continues 
to accumulate. 

With this additional procedure, unloading on any scale may be 
incorporated into the analysis. Experience has shown that the iteration 
rarely requires more than one cycle for internal consistency to be 
achieved. We may also remark that the present computer programs 
exclude Bauschinger effects largely as a matter of convenience. Had 
we the requisite data, we should be able to program not only this 
effect but cyclic stress-strain behavior as well. We recognize, 
however, that obtaining such data is itself a difficult task. 

Automatic Loading: The procedure outlined above works well 

with smooth stress-strain curves, i.e., those with a gradually changing 
slope all the way from initial yield to high strain levels. In other 



57 


situations such as that of a bilinear curve, additional steps are 
needed. The reason, of course, is that there is every chance for 
overshooting the curve at its knee when load progression is externally 
specified. 

We thus introduce the option of what is termed automatic loading. 
The analyst merely prescribes unit load steps, and the actual load 
increment level is internally scaled, once the solution is beyond the 
stage of candidacy. The scheme, operative for displacement, force, 
or mixed loading in any order, works as follows. 

Once the candidate solution, based on unit loading, is found 
acceptable, two potential scale factors are determined. One is found 
so that just one material point yields (or re-yields), and the other 
so that the increment in total octahedral strain nowhere exceeds a 
small fraction of its current value. The lesser of these two factors 
is selected and all incremental quantities are scaled by its value. 

In this manner, the stress-strain curve is tracked with no little 
precision. It may be noted, however, that the cost of solving a given 
problem is greater using automatically rather than externally determined 
load step size, since many increments are very small compared to what 
the analyst might otherwise have set. 

This method is especially useful in multi-cycle excitation in 
which progression of one system of loads is interrupted and another 
takes over in such a fashion as to force local unloading and perhaps 
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reloading. That the procedure is admissible is a direct result of the 
quasilinear nature of our original formulation; the programs in hand 
merely exploit this feature. 

An Alternate Approach: The foregoing discussion has all been 

in the context of converting differential equations to algebraic 
increment equations such as (48). While this time-honored progression 
is certainly viable, one might inquire whether an alternate procedure 
might be used whereby a functional is to be extremized. Again, the 
quasilinear nature of the formulation proves useful, for we may 
utilize such an approach in conjunction with any of the four integral 
theorems outlined above. Indeed, Jones [4] has shown that, with a 
suitable choice of minimization method, such an approach proves 
efficacious in solving elasto-plastic problems. Incorporation of 
the procedures described above, e.g., checking for local unloading, 
gives a method for problem solving that is especially useful in 
treating cases that Involve plates and shells which perforce include 
both an elastic core and a variable position of the neutral surface. 
Much of the advantage that is realized follows from the fact that 
the results of the kth increment are an excellent estimate for the 
behavior in the (k+l)th increment, and convergence of the minimization 
algorithm is quite rapid. 
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Remark on Infinite Domains: From a research standpoint, one of 

the limitations on numerical methods for solving problems is that the 
analyst is constrained to treat problems in a finite domain. It is 
almost truistic, however, that analytical results against which he or 
she would compare numerical data are for infinite domains. As a 
result, such comparisons are not always on a equal basis and 
evaluation of numerical findings is sometimes impeded. 

Riccardella has recently suggested a means for overcoming this 
limitation. Recapitulating (19)-(21) we have the (incremental) 

stress-strain relations and differential equations in the form 


6a ij /2w ' E ijU 6u k,£ 
(E ijki Su M>’j - 0 

Subject to boundary conditions on and/or 

5t i = 2yE ijk£ 6u k,£ n j 


(19) 

( 20 ) 


. (21a) 


Let us suppose now that the problem at hand involves an infinite 
domain with some sort of localized feature such as a cut-out or 
perforation. We then decompose the original problem into two parts. 

The first is the infinite domain with no local feature but subject to 
remote loading. The solution to such a problem is normally described 
in a simple manner, and we denote its solution in terms of the quantities 
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Su. and <Sa„ . In many cases of interest, the first problem is fully 
elastic and may thereby be represented analytically. 

The second problem is the residual problem for which (20) becomes 

or, if 

6uT = 6u.-Su. 

l li 


then 



(49) 


In (49) , it is important to bear in mind that depends in a 

complicated and nonlinear fashion on the stress field for the 
full problem; E^ k £ is not therefore subject to simple decomposition. 
Recognizing, however, that there is a value of E^ k £ pertinent to the 
first problem, we may rewrite (19) as 

(6a ij + 6a i:j )/2y = (E ijU + E ijk£ ) 6u^ £ + E_ k£ fiu^ 

oo 

or, since we define E^ k £ by 


5a^/2y = E. 




we have 


‘ E ijU Su k,i * E ijU 6 \,£ 


It follows that (21) leads to 

r ^ r “ 


5t i = 2u [ E ijk£ 5u k ,t * E ijk£ S \,£ ln j 


(50) 
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The second problem is thus described by the differential equation 

r r 

(49) subject to local boundary conditions such that 6u^ and/or 6t.' 

00 CO 

annihilate the residuals of 6u^ and/or <5t^. The second problem 
nominally involves an infinite domain; inasmuch, however, as the 
residual stresses and displacements die off quickly with distance 
from the local feature, a large but finite domain may be adequate 
for analysis . ■ 

We have not in fact used this procedure to date. It may be 
observed, however, that the formulation is altogether straightforward 
and follows that often used in elasticity for problems of this type, 
i.e., where there is a local disturbance in an otherwise uniform 
field. The one complexity for elasto-plastic flow — as opposed to 
elasticity — is the manner in which the generalized modulus E^^ 
is decomposed. 

Results: A variety of results has been obtained with this 

approach, most dealing with planar crack problems. Early work in 
plane stress [23] and subsequent effort in plane strain [18] was 
directed toward overall characterization of elasto-plastic response 
of cracked bodies. Concern since has focussed on deformed crack 
shapes in stretching [15,19] and bending [4]. Considerable effort 
has been directed toward the question of physical fidelity of computed 
results, and the finding is in the main satisfactory [16,17]. There 
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is significant sensitivity to the accuracy with which material data 
are represented [10], and some work has been directed toward reconciling 
experimental measurements with material response [13]. This issue is 
far from resolved and, indeed, remains as one of the primary technological 
questions of mechanics. 
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CONCLUDING REMARKS 

We have sought in this report to provide an overview of the theory 
of elasto-plastic flow without becoming deeply involved in the great 
array of detail that perforce accompanies its implementation and use. 
Our emphasis has been largely on three main aspects which deserve 
restatement at this point. First, the theory has an extensive scope. 

It can accommodate a broad variety of flow rules which reflect 
work -hardening, plastic compressibility, anisotropy, small or large 
deformations, and so on; the choice of all these resides with the 
analyst. Next, the theory is mathematically complete and coherent, 
and the analyst need have no concern for, say, the conditions under 
which it will 'work' or not. Finally, the theory is straightforward 
to use and gives results that are physically realistic insofar as we' 
have been able to judge. 

It should also be mentioned that the theory is economical to use. 
The codes now in hand, although refined over a period of some years, 
consume only modest amounts of computer storage and time, as a few 
examples will show. Using a Univac 1108, our planar codes for finite 
elements will accommodate up to 600 degrees of freedom but need less 
than 45,000 words of storage. Running times for 242 degrees of freedom 
are about 4 seconds per load step and for 436 degrees of freedom, under 
13 seconds per load step. Costs for up to 10 2 load steps are thereby 
within reach for many purposes. Indeed we have on occasion run problems 
consisting of 35 to 40 steps for classroom use. 
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The approach then is attractive in a variety of respects, and we 
look forward to its use by more people. At the same time, there are 
several items worthy of attention that should be mentioned. We 
believe that solution methods other than simple finite elements and the 
boundary-integral technique should be explored as potential vehicles 
for solving problems of elasto-plastic flow. We have looked at finite 
differences, but the context of that effort made statement of boundary 
conditions awkward. Yet there is reason to suppose that some classes 
of problem would-be tractable in finite differences. More complex 
finite elements, including special elements as might be used to 
examine locally singular behavior, merit study and possible development. 

We suggest that the breadth of the theory is in need of exploration. 
Such features as anisotropy and plastic compressibility are natural 
extensions of work done to date. Any effort along these lines, however, 
should be coupled with careful experimentation so that extraction of 
material properties procedes in parallel with development of a 
capability for solving problems. 

In addition to the features noted above, we believe that the 
economies attainable with this approach make study of cyclic loading 
especially attractive. Two important technological processes, fatigue 
response of structures and the Bauschinger effect, are thereby accessible 
to study, and we look for renewed analytical efforts along these lines. 
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The theory is also amenable to purely analytical treatments such 
as the important results reported recently by Budiansky [24]. The 
availability of solutions to non-trival problems provides both insight 
into the process of plastic flow and 'check problems' for developers 
of computer programs. Such work is usually inexpensive to perform 
but has far-reaching value. 

Finally, we urge the technical community to use this method. 

Even without further develppment and extension, it offers considerable 
advantage at reasonable cost and it provides information that may in 
some circumstances prove valuable to have. Thus incorporation of the 
theory of elasto-plastic flow into more general structural programs 
is warranted and we urge such effort upon the technical community. 
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